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EVALUATION 


This  effort  is  a  part  of  Research  Area  7,  Electronics,  Sub-Area  3, 
Communications.  The  objective  is  to  investigate  the  convergence 
properties  of  adaptive  algorithms  for  arrays  with  a  large  number  of 
antenna  elements  or  weights.  This  research  work  supports  RADC  TPO 
Sub-Thrust  4B1 ,  Adaptive  Processing  for  Communications.  The  overall 
objective  of  this  sub-thrust  is  to  advance  the  state-of-the-art  in 
adaptive  spatial  processing  to  provide  an  Electronic  Counter  Counter¬ 
measure  (ECCM)  capability  for  Air  Force  Communications  Systems. 

In  this  research  effort,  it  war.  demonstrated  that  the  judicious 
placement  of  a  few  extra  antenna  elements  provides  significant 
performance  improvement.  It  was  also  shown  that  convergence  rates  and 
beam  patterns  are  closely  related  to  the  eigenvalues  of  the  antenna/ 
signal  system.  A  closed  form  expression  was  derived  for  the  eigenvalues 
of  a  general  antenna  array  with  two  jammers.  From  this  expression,  the 
estimated  convergence  rate  and  the  optimum  beam  pattern  are  found  for 
particular  antenna  array  geometries. 

This  is  a  significant  step  in  providing  the  analytical  expressions 
enabling  selecting  the  optimum  antenna  array  geometry  for  convergence 
and  beam  pattern  shape  given  bounds  on  overall  physical  array  dimensions 
and  numbers  of  antenna  elements. 

'  '■  ,  / 

<>(  <Ll  u't  /l  u 

JOHN  A.  GRANIERO 
Project  Engineer 


I.  INTRODUCTION 


Adaptive  antenna  array  systems  can  be  used  to  steer  nulls  in  the 
direction  of  arrival  of  interfering  signal  sources  (jammers) J  Previous 
analyses  and  simulations  relating  to  these  adaptive  null  steerers  have 
focused  attention  on  the  nulling  of  a  single  jammer  or  of  multiple 
jammers  with  equal  power  and  wide  angular  separation. 

When  two  or  more  jammers  with  widely  differing  power  levels  are 
involved,  the  time  required  to  null  all  jammers  can  become  very  long  in 
comparison  with  the  time  required  to  null  a  single  jammer.  This  problem 
has  been  pointed  out  by  White®  and  Reed  zJt  at*  Furthermore,  the 
antenna  pattern  at  convergence  may  fail  to  exhibit  the  sharp  nulls  that 
are  needed  to  eliminate  interference  with  minimal  Impact  on  desired 
signals. 

The  purpose  of  the  research  reported  in  this  paper  Is  twofold. 
First,  the  impact  of  surplus  elements  on  convergence  rate  is  explored. 
Analyses  and  computer  simulations  are  presented  to  demonstrate  the 
effectiveness  of  adding  extra  antenna  elements  to  the  adaptive  null 
steerer  (ANS)  in  order  to  speed  up  the  adaptive  nulling  process  for  two 
closely  spaced  jammers  with  disparate  power  levels.  Second,  the  effect 
of  surplus  elements  on  the  converged  beam  pattern  is  considered.  A 
closed  form  expression  for  the  optimal  beam  pattern  is  used  to  directly 
compute  the  converged  beam  pattern.  This  closed  form  expression 
provides  confirmation  of  the  simulation  results  and  furnishes  a  computa¬ 
tionally  efficient  means  of  exploring  the  converged  beam  patterns  for 
different  element  configurations. 


II.  DEMONSTRATION  OF  THE  PROBLEM 


A.  Rate  of  Convergence 

To  demonstrate  the  slowness  In  the  nulling  process,  a  minimal 
adaptive  sldelobe  canceller  structure  (Fig.  2-1)  having  one  primary 
element  and  two  auxiliary  elements  Is  employed  to  null  two  jammers* 

designated  as  follows: 

J.  *  *  Jammer  1  with  power  100,  Incident 
1  at  40°,  frequency  *  0.231/T 

J9  =  Jammer  2  with  power  1,  Incident 
c  at  55°,  frequency  *  0.170/T 

T  »  sampling  period 

The  array  elements  are  spaced  a  quarter  wavelength  apart  at  the  frequency 
of  0.2/T.  The  angular  position  of  J^  Is  chosen  rather  arbitrarily,  while 
that  of  J2  Is  chosen  reasonably  close  to  Jj.  The  desired  signal*  not 
explicitly  included  In  this  simulation,  was  assumed  to  be  of  much  lower 
power  level  than  that  of  either  jammer.  No  constraints  were  used  in  this 
system  to  prevent  cancellation  of  the  desired  signal.  The  intent  is  to 
obtain  a  beam  pattern  of  unit  gain  In  all  directions  except  where  high 
power  jammers  are  Incident.  There,  the  gain  should  be  very  small  indeed. 

To  focus  attention  on  the  essence  of  the  problem  It  is  assuaed 
that  the  jammers  are  narrowband  signals  which  may  be  approximated  by 
sinusoids.  Thus  90°  phase  shifters  (Instead  of  tapped  delay  lines)  are 
used  on  each  auxiliary  antenna  element.  The  LMS  adaptive  algorithm 
Is  used  In  a  simulation  Involving  the  configuration  of  Fig.  2-1.  After 
a  few  adaptations  the  weights  are  frozen  and  the  resulting  antenna  beam 
pattern*  Is  plotted.  Then  the  adaptive  process  Is  resumed  and  the 

*The  beam  pattern  Is  the  array  response  to  a  unit  power  test  signal 
which  is  swept  over  360°  while  keeping  the  adjustable  weights  frozen. 
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Auxiliary  Elements 


pattern  plotting  procedure  is  repeated  at  frequent  intervals.  It  may 
be  noted  that  no  noise  is  added  to  the  signals  in  order  to  concentrate 
exclusively  on  the  convergence  problems  due  to  the  jamming  signals. 

Figures  2-2a  through  2-2d  show  a  time  sequence  of  beam  patterns 
achieved  by  the  above  simulation  as  adaptation  proceeds.  The  arrows  In 
the* figure  indicate  the  positions  of  the  jammers. 

From  Fig.  2-2a  it  is  clear  that  the  high-power  jammer  is  nulled 
quickly,  but  the  low-power  jammer  is  essentially  unaffected  at  the 
outset.  In  Fig.  2-2b  the  pattern  has  changed  substantially;  but  it  still 
shows  little  nulling  effect  in  the  area  of  the  low-power  jammer.  After 
about  40,000  adapts  the  low-power  signal  begins  to  be  nulled  as  shown 
in  Fig.  2-2c.  Well  formed  nulls  begin  to  occur  after  about  70,000 
adapts.  Good  nulls  are  achieved  by  about  120,000  adapts.  For  the 
purpose  of  this  paper,  convergence  may  be  defined  as  the  minimum  number 
of  adapts  needed  to  reduce  the  output  power  to  below  1  x  10“®.  For  the 
configuration  in  Fig.  2-1,  this  occurs  around  141,000  iterations.  The 
converged  pattern  is  shown  in  Fig.  2-2d.  If  the  sampling  rate  were 
1  MHz,  141,000  adapts  would  require  0.141  sec. 

The  sluggishness  of  the  adaptive  process  in  nulling  the  low-power 
jammer  may  be  explained  with  the  aid  of  the  analysis  presented  in 
section  III. 

B.  Final  Pattern  Characteristics 

The  converged  beam  pattern  (Fig.  2-2d)  of  the  minimal  system  is 
symmetric  about  0°  because  a  line  array  of  antenna  elements  is  unable 
to  distinguish  between  positive  and  negative  angles  of  arrival.  The 
additional  nulls  introduced  by  symmetry  are  undesired  since  they  may  be 
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NOTE:  Beam  patterns  plotted  for  f  ■  .2/T 


Fig.  2-2.  A  time  sequence  of  beam  patterns 
for  the  ANS  of  Fig.  2-1. 


In  the  direction  of  the  desired  signal.  The  broadness  of  the  nulls 
also  may  cause  attenuation  of  any  signal  with  an  angular  location  that 
is  close  to  that  of  either  of  the  symmetric  nulls.  A  desirable  pattern 
would  be  omnidirectional  with  sharp  nulls  in  only  the  directions  of  the 
jammers . 

The  improvements  in  the  final  pattern  due  to  restructuring  of  the 
antenna  array  will  be  detailed  in  section  VII. 


III.  ANALYSIS  RELATING  TO  CONVERGENCE  RATE  OF  ADAPTIVE  ARRAYS 


A  generalized  version  of  the  adaptive  antenna  array  of  Fig.  2-1  is 
shown  in  Fig.  3-1.  The  kth  auxiliary  element  is  located  at  a  distance 
t^  (given  in  units  of  wavelengths)  from  the  primary  element  at  aa  angle 
with  an  arbitrary  reference  axis.  Two  jammers  (J^  and  Jg)  with 
powers  p-j  and  P2  (at  each  antenna  element)  arrive  at  angles  ^  and  ^ 
as  shown  in  the  figure. 

Central  to  the  analysis  of  the  adaptive  linear  combiner  of  the  type 
used  in  Fig.  3-1  is  the  correlation  matrix  R  (Widrow  &£  at )7  defined  by 

R  -  E(X1  ,Xj)  (3-1) 

where  E(X^  )  is  the  statistical  correlation  between  the  signals 

X..  and  Xj  where  i ,  j  =  1 ,  2 . 2n.  The  speed  of  convergence  is 

intimately  related  to  the  eigenvalues  of  R,  as  can  be  seen  from  the 
following. 

For  stability  of  the  LMS  adaptive  algorithm  [see  (24)  of  Widrow7) 
the  adaptation  constant  y  must  satisfy 

IMmax  >u>0  (3-2a) 

where  X  is  defined  as  the  maximum  eigenvalue  of  the  R  matrix, 
max 

Maximum  speed  is  obtained  with 

P  =  l/2Xmax  .  (3-2b) 

The  dynamics  of  the  adaptive  system  depend  only  on  the  nonzero 
eigenvalues  of  R.  For  every  nonzero  eigenvalue  Xj,  the  time  constant 
t,j  of  the  corresponding  mode  is  given  by  [see  (27)  of  Widrow7) 
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(3-3) 


Tj  =  l/2yXi 

(t.  Is  dimensionless.  It  is  expressed  in  number  of  adapts.) 

The  rate  of  convergence  with  a  specified  pis  thus  bounded  above  by 
the  minimum  nonzero  eigenvalue  A^.  From  (3-3): 

tmx  ■  l/2l,X.1ii  •  (S-4) 

Using  relation  (3-2b),  when  adapting  at  the  maximum  rate  consistent  with 
the  requirements  for  stability,  the  time  constant  of  the  slowest  mode 
will  be 

Tmax  ~  Nnax^Nnin  *  (3-5) 

In  practice  the  p  Is  required  to  be  small  compared  with  1/A  to  keep 
the  misadjustment+  within  an  acceptable  limit.  Therefore  is 
greater  than  that  Indicated  by  equation  (3-5).  Nevertheless,  the  ratio 

\iax^Am1n  *n  the  r*9ht  hand  side  (RHS)  (3-5)  Is  a  useful  measure  for 
comparing  speeds  of  convergence,  as  will  be  seen  from  the  examples 
discussed  In  sections  IV  and  V. 

From  a  detailed  analysis  given  In  Appendix  A,  It  is  shown  that  the 
eigenvalues  for  the  general  case  in  Fig.  3-1  take  on  the  values  given 
by  the  following  expressions: 


+Us1ng  the  definition  from  Reference  7: 

A  M  Average  excess  mean  squared  error  (mse) 

Mlsadjustment  H  mse  obtained  with  optimum  weight  setting 


V|  =  2'(P'j+p2)  +  (pi"p2^  +  P2P1 1 ^*2 1 

^2  =  ?^Pl+p2^  ~  ^7  ^pl”p2^  +  P2P1 1 ^1  ^2 ^ 


Xn+1  =  Xl;  Xn+2  =  X2 

*  0  for  m  t  {1 ,2,n+l ,n+2) 

where 


j  -  1,2  (jammer  Index) 

U*  =  Complex  conjugate  transpose  of  Uj 


♦«  5  2"\  cos<*j",k) 

k  *  1,  ...»  n  (auxiliary  element  index) 

t .  =  distance  of  kth  antenna  element  from  the 
*  primary  antenna  element  (ty  given  In  units 

of  wave  length) 


a.  =  Angle  of  kth  element  with  respect  to  the 
*  reference  axis  (see  Fig.  3). 

^  »  Angle  of  1th  jammer 
p.|  ■  Power  of  1th  jammer 


- 1 
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The  ratio  of  maximum  eigenvalue  to  minimum  nonzero  eigenvalue  is 


thus  given  by 


max  _  1 

\nin  X2 


(P1+P2) 


(prp2)‘ 


+  plp2 


.  *  .2 


n  J 


(pl +p2^ 


(prp2)Z ,  P,p2  l«iV 


(3-7) 


From  the  definition  of  in  (3-6)  it  Is  seen  that  [U*U2[  varies 

between  0  and  n.  Thus  for  given  and  p2  (p-j  >  p2),  the  ratio 

Am  /A  varies  between  p./p0  and  «.  Am  /A ,  .  must  be  kept  small  in 
max  ip  in  l  c.  max  min 

order  to  keep  t  small  [see  (3-5)].  Therefore  the  number  of  elements 

IMaX 

and  their  geometrical  arrangement  in  the  array  are  to  be  chosen  such 
that  the  ratio  ^^^min  is  kept  near  P^/q2  for  most  and  ^2  of 
practical  interest.  Note  that  ^  ^2  ^  Uj  U2  => 

|U*U2|  -►  n, which  implies  from  (3-7)  that  Amax/Am^n  ®. 

To  obtain  a  better  understanding  of  the  implications  of  (3-6)  and 
(3-7),  compare  two  special  classes,  namely: 

(i)  P  •  P,  -  P2 

Substituting  in  (3-6)  and  (3-7)  yields 

xi  ■  Vi  ■  p<n  *  Ku*i> 

x2  *  xn+2  ’  p("  "  I U*U2 1  * 

x  •  0  otherwise 
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(3-9) 


jnax  =  1  + 

Xm1n  1  -  |U*U2|/n 


(11)  p1  »  p2 


Making  first-order  approximations  in  (3-6)  and  (3-7)  it  cam  be 
seen  that 


Xn+1  =  npl 


Xn+2  *  np2 


[■¥ 


\  *  0  otherwise 

m 


(3-10) 


max 

^in 


U— 3, - 1  •  [ - 1 - 1  .  (3-11) 

2  [1  +  |U|U2|/n J  |_1  -  lU^I/nJ 


When  |U*U2|/n  approaches  1  It  is  seen  from  (3-9)  and  (3-11)  that 
the  eigenvalue  disparity  (i.e.,.the  ratio  \,ax/\„fn)  Increases  without 
bound,  regardless  of  the  disparity  in  jammer  power  levels  (l.e.,  the 
ratio  p-j/p2).  The  eigenvalue  disparity  in  (3-11)  is  approximately 
0.25  P|/p2  times  the  eigenvalue  disparity  in  (3-9).  It  may  be  observed 
from  (3-10)  that  the  high-power  jammer  determines  the  largest  eigenvalue. 
This  implies  [according  to  (3-3)]  that  the  high-power  jammer  Is  attenuat¬ 
ed  quickly.  From  (3-10)  it  is  also  seen  that  the  smaller  eigenvalue  Is 
determined  by  the  low-power  jammer,  by  the  angular  separation  between 
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the  jammers,  and  by  the  array  configuration.  Thus  for  (U^l/n 

approaching  1,  we  see  that  this  eigenvalue  becomes  quite  small.  This 

explains  why  the  low-power  jammer  was  attenuated  very  slowly  in  Fig.  2-2 
* 

When  I U-|U2 1 /n  approaches  0  it  is  seen  from  (3-8)  that  both  the 

eigenvalues  are  the  same  and  hence  the  ratio  A  .  approaches  1, 

max  min 

the  same  ratio  in  class  (ii)  becomes  only  p^/pg. 

To  summarize  the  above,  it  may  be  stated  that  the  eigenvalue 

* 

disparity  varies  between  p^/p2  and  «  as  the  |U^U2 [ /n  varies  between  0 
and  1.  | U-j l>2 1 /n  in  turn  depends  upon  the  angular  positions 

the  two  jammers  and  on  the  number  of  antenna  elements  and  their 
geometrical  arrangement.  For  a  given  range  of  values  of  and  ^ 
the  number  of  elements  and  their  arrangement  may  be  chosen  such  that 

i*i 

| ll-j U2 1 /n  may  be  kept  small,  which  in  turn  will  keep  the  eigenvalue 
disparity  near  the  jammer  power  disparity. 

Examples  of  various  array  configurations  are  examined  in  the  next 
section  in  order  to  study  their  effect  on  the  eigenvalue  disparity. 


IV.  ESTIMATING  RELATIVE  CONVERGENCE  RATES  IN  VARIOUS  ARRAY  CONFIGURATIONS 

Several  antenna  array  configurations  were  tried  in  order  to  study 
the  eigenvalue  disparity  for  different  numbers  of  elements  and  various 
element  configurations.  One  of  these  configurations  has  already  been 
shown  in  Fig.  3-1.  where  the  element  locations  are  shown  together  with 
the  phase  shifters  and  the  adaptive  linear  combiner.  Other  array  con¬ 
figurations  are  shown  in  the  "a"  parts  of  the  figures  in  this  section. 

In  these  latter  figures,  only  the  element  locations  are  shown  explicitly; 
phase  shifting  and  adaptive  linear  combining  are  understood  to  follow  the 
pattern  established  in  Figures  2-1  and  3-1. 

Figure  4-1  and  the  “b"  parts  of  the  remaining  figures  in  the  section 
are  plots  of  the  normalized  eigenvalue  disparity  (NED),  which  is  defined 

as  ^\iax^\nin^pl^p2^  for  the  Partlcu^ar  configuration  under  considera¬ 
tion.  In  each  case  the  angular  position  of  jammer  1,  i.e.,  is  fixed 
arbitrarily  at  40°  while  that  of  jamner  2,  I.e.,  is  varied  between 
-180°  and  180°.  Normalized  eigenvalue  disparity  is  plotted  as  a  function 
of  ip2. 

Figure  4-1  is  the  NED  plot  for  the  linear  array  with  two  auxiliary 
elements  that  was  used  in  Section  II  to  show  that  very  slow  convergence 
can  be  experienced  for  closely  spaced  jammers  with  disparate  power 
levels.  The  particular  case  that  was  examined  involved  a  55°  incidence 
angle  for  the  second  jammer  Jg.  From  Figure  4-1  It  can  be  seen  that  the 
NED  for  this  Incidence  angle  Is  approximately  45.  This  large  eigenvalue 
disparity  shows  why  a  very  long  adaptation  period  was  required  for  this 
particular  antenna/signal  configuration. 


Figure  4-1  also  indicates  that  similar  convergence  problems  can  be 


anticipated  over  a  broad  range  of  possible  ^  arrival  angles. 

Particularly  difficult  convergence  problems  are  indicated  at  arrival 
angles  that  approximate  the  arrival  angle.  This  result  runs  contrary 
to  Intuition  and  experience  since  it  would  be  expected  that  a  single, 
slightly  broadened  null  would  form  quickly  and  eliminate  both  jammers  In 
such  cases.  Appendix  D  examines  the  apparent  conflict  between  the 
results  indicated  by  the  NED  plot  and  the  results  that  have  been  observed 
in  practice.  It  is  shown  that,  almost  paradoxically,  both  viewpoints 
have  their  validity  and  that  the  choice  of  viewpoint  depends  upon  the 
specification  of  initial  conditions  and  the  definition  of  convergence. 

The  NED  plot  for  the  linear  configuration  of  Fig.  4-2a  is  shown  in 
Fig.  4-2b.  This  plot  is  very  similar  to  that  of  Fig.  4-1,  though  the 
value  of  NED  for  each  i|u  is  smaller  than  the  corresponding  value  in 
Fig.  4-1.  At  =  55°,  NED  for  the  three  auxiliary  case  is  17.1,  which 
is  2.6  times  smaller  than  the  NED  in  the  linear  configuration  with  two 
auxiliaries.  This  Indicates  that  the  nulling  process  in  the  array  with 
three  auxiliaries  will  be  correspondingly  faster  than  that  in  the  array 
with  two  auxiliaries. 

A  simple  way  to  break  the  symmetry  of  the  linear  array  is  to  use  a 
triangular  array  such  as  that  shown  in  Figure  4-3a.  Comparing  Figure 
4-3b  with  Figs.  4-1  and  4-2b,  it  is  evident  that,  in  the  region  where 
small,  the  normalized  eigenvalue  disparity  (NED)  is  sub¬ 
stantially  less  for  the  triangular  configuration  than  the  NED  for  the 
linear  configurations  with  two  or  three  auxiliary  elements.  As  a 
specific  example,  at  i if>2  =  55°  the  NED  for  the  triangular  configuration 
is  12.5,  while  the  NED  for  the  linear  configuration  (Fig.  4-1)  is  45.0. 
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Figure  4-2.  Antenna  Geometry  and  NED  Plot  for  the 

Linear  Array  with  Three  Auxiliary  Elements. 
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Figure  4-3.  Antenna  Geometry  and  NEO  Plot  for  the  Triangular  Array. 
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That  is,  the  NED  in  the  triangular*  configuration  is  3.6  times  smaller 
than  the  NED  in  the  linear  configuration.  Using  the  lower  bound  on 
x  [i.e.,  (3-5)]  it  is  estimated  that  the  adaptive  nulling  process  in 
the  triangular  configuration  is  about  3.6  times  faster  than  that  in  the 
linear  configuration. 

The  hexagonal  configuration  of  Fig.  4-4a,  a  natural  extension  of 
the  triangular  configuration  just  described,  was  also  tried.  In  the 
region  -50°  <  \i>  <  110°,  the  NED  plot  shown  in  Fig.  4-4b  appears  identi¬ 
cal  to  that  in  Fig.  4-3b.  In  the  other  regions  the  NED  is  somewhat 
better  (i.e.,  smaller)  than  that  in  Fig.  4-3b.  Thus  the  hexagonal  con¬ 
figuration  may  be  expected  to  exhibit  convergence  behavior  much  like  the 
triangular.  At  ^  -  55°  the  NED  is  the  same  as  in  the  triangular  case 
and  hence  the  convergence  rate  in  both  cases  is  expected  to  be  the  same. 

Next  consider  the  configuration  shown  in  the  Fig.  4-5a,  which 
consists  of  two  concentric  triangles.  The  inner  triangle  is  the  same 
as  that  of  Fig.  4-3a.  The  outer  triangle  has  its  vertices  at  a  radius 
of  7.75  wavelengths  (which  is  30  times  the  radius  of  the  inner  one). 

The  resulting  NED  plot  is  shown  in  Fig.  4-5b.  Except  for  a  very  small 
region  around  40°  (the  position  of  Jammer  1)  the  value  of  NED  rarely 
exceeds  10.  For  the  example  of  ^  =  55°,  the  value  of  NED  is  1.13. 

This  values  is  approximately  11  times  smaller  than  the  corresponding  NED 
for  the  single  small  triangle  case  (Fig.  4-3b).  This  indicates  a  corre¬ 
spondingly  faster  convergence  rate  in  the  two  triangle  case. 

Comparing  the  two-triangle  case  (Fig.  4-5b)  with  the  two-auxiliary- 
element  linear  configuration  (Fig.  4-1)  the  NED  in  the  two- triangle  case 
Is  40  times  smaller  than  the  corresponding  NED  in  the  linear  case.  Thus 
it  appears  that  the  nulling  process  in  the  two  triangle  case  will  be  40 
times  faster  than  that  in  the  linear  case. 
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The  reasoning  which  led  to  the  selection  of  the  two-triangle  case 
for  study  was  as  follows:  It  may  be  seen  that  elements  spaced  at  a  long 
distance  from  the  primary  element  experience  a  large  time  difference  in 
the  signal.  Hence  intuitiveZy  a  lajige  separation  between  eZements 
(a  large  aperatune)  may  tend  to  magnify  the  apparent  angular  separation 
between  closely  spaced  Jammers.  This  magnification  tends  to  reduce  the 
disparity  between  eigenvalues  of  the  adaptive  null  steerer. 

As  a  final  case  for  examining  relative  convergence  rates,  the  three- 
triangle  configuration  of  Fig.  4-6a  was  considered.  This  is  an  extension 
of  the  two-triangle  case;  the  inner  and  outer  triangles  of  elements 
remain  unchanged.  Figure  4-6b  is  the  NED  plot  for  this  configuration. 

The  NED  values  in  this  case  fall  generally  below  the  NED  values  for  the 
two-triangle  configuration,  and  the  strong  spike  at  about  25°  has  been 
suppressed.  The  improvement  is,  however,  relatively  small  at  most 
angles.  In  particular,  the  NED  value  at  ^  51  55°  *s  an  insignifi¬ 

cant  improvement  over  the  value  of  1.13  that  was  obtained  with  the  two- 
triangle  arrangement.  It  appears  that  convergence  for  the  three- 
triangle  case  will  be  about  41  times  faster  than  the  convergence 
performance  exhibited  by  the  linear  array  with  two  auxiliaries. 

The  eigenvalue  disparities  have  been  discussed  for  various  array 

configurations  based  upon  arrangements  of  elements  in  simple  linear  and 

triangular  patterns.  The  eigenvalue  disparity  Is  related  to  the  time 

constant  t  „  of  the  slowest  mode  when  adapting  with  the  value  of  p 
max 

giving  the  highest  convergence  rate,  in  accord  with  equation  (3-5).  The 
results  obtained  from  the  simulations  given  In  the  following  section 
confirm  that  the  eigenvalue  disparity. can  provide  useful  Information 
about  convergence  rate. 
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V.  SIMULATION  RESULTS 


The  principal  thrust  of  this  section  will  be  the  discussion  of 
convergence  rate  and  a  comparison  of  simulation  results  with  the 
estimates  from  the  foregoing  section.  The  beam  patterns  at  convergence 
will  also  be  documented  in  this  section  for  each  antenna  configuration. 
Discussion  of  the  beam  patterns  will,  however,  be  deferred  until 
Section  VII  where  the  discussion  can  Include  material  related  to  the 
optimal  beam  pattern  as  well  as  the  simulation  results. 

Simulations  for  the  cases  discussed  in  Section  IV  were  conducted 
in  a  manner  similar  to  that  described  in  Section  II  for  the  linear  array 
with  two  auxiliary  elements.  The  same  simulation  program  was  applied 
to  all  the  arrays,  and  the  number  of  adapts  required  to  reduce  output 
power  below  1  x  10~6  was  determined  for  each  array  configuration. 
Additionally,  the  beam  pattern  at  convergence  was  plotted  for  each  case. 

Much  of  the  data  concerning  convergence  rate  characteristics  for 
the  different  array  configurations  is  summarized  in  Table  1.  The 
columns  at  the  left  side  of  the  table  describe  the  array  configurations 
and  give  cross  references  to  the  figures  that  illustrate  each  geometry. 
The  center  section  of  the  table  brings  forward  the  pertinent  information 
from  the  normalized  eigenvalue  disparity  (NED)  curves  of  Section  IV  for 
the  particular  case  of  interest,  i.e.,  for  jammer  azimuths  of  40°  and 
55°.  Additionally,  this  section  shows  the  relative  convergence  rates 
with  the  baseline  configuration  taken  as  the  linear  array  with  two 
auxiliary  elements.  The  right  side  of  the  table  gives  a  summary  of  the 
simulation  data.  The  adaptation  constant  is  shown  for  each  case 
together  with  the  observed  number  of  adaptations  for  convergence  and 
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CONVERGENCE  RATE  CHARACTERISTICS  OF  VARIOUS  ARRAY  CONFIGURATIONS 
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an  observed  relative  convergence  rate.  Once  again,  the  relative 
convergence  rate  Is  based  upon  the  convergence  rate  of  the  linear  array 
with  two  auxiliary  elements. 

Appropriate  selection  of  the  adaptation  constant  y  was  an  Important 
consideration  In  conducting  the  simulation  experiments.  The  major 
concern  was  to  assure  that  convergence  times  could  be  meaningfully 
compared  for  the  various  array  configurations;  this  was  accomplished  by 
setting  the  adaptation  constant  experimentally  for  the  linear  array 
with  two  auxiliary  elements,  and  then  adjusting  the  adaptation  constant 
for  other  cases  so  that  the  product  of  y  and  the  array  Input  power 
remained  a  constant.  The  value  of  y  *  1  x  10  was  selected  for  the 
linear  array  with  two  auxiliary  elements  by  conducting  a  few  experimental 
runs.  A  smaller  value  of  y  resulted  In  a  slower  convergence  rate  than 
that  shown  earlier  In  Figure  2-2,  with  no  Improvement  In  the  pattern.  A 
larger  value  of  y  resulted  in  poor  beam  patterns  due  to  Increased  mls- 
adjustment  In  the  weight  vector  [see  Hldrow  8].  For  other  array  con¬ 
figurations  listed  In  the  table,  y  was  set  Inversely  proportional  to  the 
number  of  auxiliary  elements  In  each  case.  If  y  were  set  smaller  than 
shown  In  the  table,  the  adaptive  nulling  process  would  be  slower  than 
that  obtained  with  the  given  y,  without  any  Improvement  in  the  team 
pattern.  If  y  were  set  larger,  the  beam  pattern  would  deteriorate  fn 
comparison  with  that  obtained  with  the  given  y. 

From  Table  1  It  Is  seen  that  the  linear  configuration  with  three 
auxiliaries  converged  2.8  times  faster  than  the  case  In  Fig.  2-1.  This 
Is  slightly  faster  than  that  estimated  on  the  basis  of  elgenvalae 
disparity.  The  converged  beam  pattern  (Fig.  5-1)  is  similar  to  that 
In  Fig.  2-2d.  The  pattern  shows  similar  notches  in  the  jammer  positions, 

l 

and  it  is  slightly  better  In  other  regions. 
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Converged  Beam  Pattern  for  the  Linear  Arrav 
with  Three  Auxiliary  Elements  (Figure  4-2a) 


For  the  triangular  array  configuration,  the  observed  convergence 
rate  Is  5.4  times  faster  than  that  for  the  linear  configuration  of 
Fig.  2-1.  It  may  be  noted  that  the  observed  convergence  rate  Is  sub¬ 
stantially  faster  than  that  estimated  In  Section  IV  by  comparison  of 
the  eigenvalue  disparity.  The  converged  beam  pattern  Is  shown  In 
Fig.  5-2. 

The  results  of  simulation  of  the  hexagonal  configuration  are 
essentially  the  same  as  those  for  the  triangular.  The  convergence  rate 
Is  the  same  as  for  the  triangular  case.  The  beam  pattern  (Fig.  5-3)  is 
strikingly  similar  to  that  shown  in  Fig.  5-2. 

The  observed  convergence  rate  in  the  two-concentric- triangle 
configuration  is  about  34  times  faster  than  that  in  the  case  in  Fig.  2-1. 
This  rate  Is  slightly  lower  than  estimated  on  the  basis  of  eigenvalue 
disparity.  The  converged  beam  pattern  Is  shown  in  Fig.  5-4. 

The  final  configuration  that  was  examined  was  the  three-concentric- 
triangle  case.  The  estimates  of  convergence  rate  Indicated  that 
performance  would  be  only  slightly  better  than  that  given  by  the  two- 
concentric-triangle  case.  In  the  simulation  a  minor  Improvement  of  100 
adaptations  was  obtained  over  the  number  of  adaptations  needed  In  the 
two-concentrlc-trlangle  case.  The  beam  pattern  for  the  configuration 
of  three  concentric  triangles  Is  shown  In  Fig.  5-5. 


ADRPTS 


Figure  5-4.  Converged  Beam  Pattern  for  the  Array  with  Auxiliary  Elements 
Configured  as  Two  Concentric  Triangles  (Figure  4-5a). 


Figure  5-5.  Converged  Beam  Pattern  for  the  Array  with  Auxiliary  Elements 
configured  as  Three  Concentric  Triangles  (Figure  4-6a). 


VI.  DERIVATION  OF  OPTIMAL  BEAM  PATTERN 

A  closed  form  solution  for  the  beam  pattern  of  the  adaptive  antenna 
array  of  Fig.  3-1  will  now  be  found. 

Let  W  be  the  weight  vector  defined  as 

W  =  [w1  w2  w3  ...  w2n]T  (6-1) 

where  w^  and  w^+n  (i=l...n)  represent  the  in-phase  and  quadrature  phase 
weights  at  element  i.  Let  X  be  the  vector  of  received  signals  at  each 
element.  Then 

X  =  [Xj  x2  ...  xR  ...  x2n]T  (6-2) 

where  the  xi  and  *(i+n)  (i ■!...«)  represent  the  in-phase  and  quadrature 
phase  components  of  the  signal  at  element  i.  Also  let  Xq  be  the  signal 
with  power  p  at  the  primary  element. 

The  output  of  the  antenna  array  is  given  as 

y  *  Xq  -  WT-X  .  (6-3) 

Consequently,  the  output  power  becomes 

E[y2]  =  p  -  2WTE[xqX]  +  WTE[X*XT3  W  .  (6-4) 

Def i  ne 

P  *  E[x0-X] 
and 

R  -  E[X.XT] 
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Derivation  of  the  sensitivity  pattern  for  the  array  depicted  in 
Figure  3-1  is  achieved  as  follows.  Assume  a  sinusoidal  plane  wave  of 
power  p  and  angle  ip  incident  on  the  array.  Then  it  can  be  shown  that 


where 


Therefore 


and 


E[Xj-xk]  =  cosUj  -  <|>k) 


(6-6) 


0  k  =  0 

2n-£kcos(ak  -  ip)  .  t  <_  k  <  n 

2Tr*£knsin(akn  -  ip)  n  +  1  £  k  £  2n 


(6-7) 
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R(^)  =  =  E 


x1x1  XjX2  .  .  .  x,x2n 


x2nxl  ‘  *  *  x2nx2n 


(6-9) 


Hence,  given  the  geometry  of  the  array  and  the  arrival  angle  of  the 
test  signal,  P  and  R  can  easily  be  calculated. 
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In  Appendix  B  it  is  shown  that 


Z,  I  Z 

k=l  Xk  °N 


jammer 


(6-10) 


where  is  the  k**1  eigenvector  and  X^  is  the  k*'*1  eigenvalue  of 

2 

Rjammer*  as  defined  ih  Appendix  C.  Recall,  also,  that  oN  is  the  power 
of  the  isotropic  white  noise  incident  upon  every  element.  The  vector 
Pjammer*  s^m^ar  structure  to  the  P  defined  in  Eq.  (6-8). 
Specifically,  it  is  determined  by  M  independent  jammers,  in  the  follow¬ 
ing  manner. 


M 


jammer 


-I 


P(^) 


i=l 


(6-11) 


where  ^  is  the  arrival  angle  of  jammer  i.  P(^)  is  defined  by 
Eqs.  (6-5, 6-6, 6-7).  Substituting  (6-10)  into  (6-3)  yields 


E[y2] 


P 


PT 

jammer 


L  ■;  ~  77 

k=l  xk  0N 


^ jammer 


(6-12) 


+  P 


jammer 


V  Qk’Qk 
2,  2 
k=l  Xk  0N 


2  Qk'°kTj 
k=l  Xk  +  °N 


jammer 


Once  the  eigenvectors  and  eigenvalues  are  established,  Eq.  (6-12) 
specifies  the  optimal  beam  pattern.  Note  that  the  above  analysis  is 
true  for  M  jammers.  Appendix  A  derives  eigenvectors  and  eigenvalues  for 
the  two  jammer  case. 
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VII.  BEAM  PATTERN  CHARACTERISTICS 


To  compare  the  converged  beam  patterns  of  different  antenna 
configurations,  a  desired  beam  pattern  must  be  established.  For  this 
study,  the  desired  beam  pattern  is  unity  gain  in  all  directions  except 
where  a  jammer  is  incident  on  the  array,  where  the  gain  should  approach 
zero. 

In  section  VI,  a  closed  form  expression  for  the  converged  beam 
pattern  was  derived.  That  expression  was  used  to  develop  a  program  for 
directly  generating  the  converged  pattern  for  a  particular  array.  The 
simulation  results  of  Fig.  2-2d,  took  141,000  adapts  to  converge.  This 
represents  a  substantial  amount  of  computer  time.  The  closed  form 
solution  was  used  to  generate  the  pattern  of  Fig.  7-la,  verifying  the 
simulation  results.  Much  less  computer  time  was  required  to  generate 
this  pattern  than  to  generate  the  simulation  results.  The  program  using 
the  closed-form  expression  allowed  quick  comparison  of  converged  beam 
patterns  for  different  array  configurations. 

Converged  beam  patterns  for  two  basic  antenna  geometries,  linear 
and  triangular,  will  be  discussed.  The  linear  array  consists  of  omni¬ 
directional  antenna  elements  placed  on  a  line.  As  mentioned  previously, 
a  linear  array  results  in  undesired  symmetric  nulls.  By  adding  additional 
equally  spaced  elements,  extra  degrees  of  freedom  are  introduced  into 
the  system.  As  can  be  seen  in  Figures  7-la,  b,  c  the  additional  array 
elements  cause  the  nulls  to  become  sharper,  and  Introduce  increasing 
ripple  into  the  beam  pattern.  The  gain  in  the  region  between  the  two 
jammers  comes  closer  to  unity  as  the  number  of  elements  is  increased. 


M  i  nwninii 
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Patterns  at  Convergence  for  Several  Linear  Arrays, 


The  triangular  arrays  consist  of  antenna  elements  positioned  as 
shown  In  Figures  4- 3a  (3  elements),  4-9a  (6  elements),  and  4-6a 
(9  elements).  Each  of  the  latter  two  figures  builds  upon  Its  predecessor 
by  adding  three  additional  elements.  Note  in  Figures  7-2a,  b,  c  that 
the  undesirable  symmetric  nulls  are  no  longer  present.  As  with  the 
linear  arrays,  as  the  number  of  elements  is  Increased,  the  nulls  sharpen, 
and  ripple  Is  introduced. 

Generally,  increasing  the  number  of  auxiliary  antenna  elements 
sharpens  the  nulls,  in  addition  to  increasing  the  convergence  rate.  A 
triangular  array  eliminates  the  unwanted  symmetric  nulls. 


Patterns  at  Convergence  for  Several  Triangular  Arrays. 


VIII.  CONCLUSIONS  AND  DIRECTIONS  FOR  FURTHER  STUDY 


It  has  been  demonstrated  that  the  addition  of  a  few  surplus 
auxiliary  elements  can  effect  a  substantial  improvement  in  the  convergence 
rate  of  an  adaptive  antenna  array  subjected  to  two  closely  spaced  jammers 
of  disparate  power  level.  A  theoretical  analysis  showed  that  the 
normalized  eigenvalue  disparity  provided  a  good  measure  of  adaptation 
rate  for  any  given  configuration.  Computer  simulations  confirmed  the 
validity  of  the  analysis  and  furnished  specific  illustrations  of  array 
geometries  with  convergence  rates  far  better  than  that  exhibited  by  a 
minimal  linear  array.  Additional  work  is  needed  to  provide  more  specific 


design  guidance  for  arrays  with  good  convergence  properties.  It  would 


be  useful  to  examine  the  question  of  optimal  convergence  rate  for  the  two- 


jammer  case,  given  some  set  of  constraints  on  array  dimensions  and  element 


count. 


Considerable  progress  has  been  made  in  developing  a  technique  for 
studying  the  impact  of  antenna  geometry  on  the  converged  beam  pattern  in 
the  two-jammer  case.  The  closed-form  expression  for  computation  of  the 
beam  pattern  provided  an  efficient  way  to  determine  the  pattern  for  any 
element  geometry  that  might  be  under  consideration.  Further  work  is 
needed,  however,  to  provide  better  Insight  into  the  synthesis  of  arrays 
with  good  beam-pattern  properties. 

The  concepts  and  results  of  this  paper  could  also  be  extended  to  the 
3-dimensional  case,  where  jammers  and  array  elements  are  no  longer 
coplanar.  This  will  require  techniques  for  representing  the  3-dimensional 
beam  pattern  In  2  dimensions.  Consideration  of  the  3-dlmensional 
extensions  has  begun. 
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Extension  of  these  results  to  arrays  containing  many  more  than 
three  or  four  elements  will  be  undertaken.  New  approximation  methods 
will  be  required  to  keep  the  algebraic  expressions  tractable  and 
insightful . 
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APPENDIX  A 


The  purpose  of  this  Appendix  Is  to  obtain  explicit  expressions  for 
the  eigenvalues  and  eigenvectors  of  the  adaptive  antenna  array  of  Fig.  3. 

Let  the  two  jammers  arriving  at  ^  and  ^  be  sinusoids.  Let  the 
weights  attached  directly  to  the  antenna  elements  be  labeled  1  through 
n,  and  those  attached  after  the  -90°  phase  shift  be  correspondingly 
labeled  n  +  1  through  2n.  See  Fig.  3. 

The  signal  at  the  Input  to  kth  weight  Is  given  by 


xk  ~  2p.j  cos(u1t+0i-«j)kl ) 

+  2p2  cos(<^t+e2-<j>k2) 


where,  for  1  *  1,  2  we  have 


Pj  *  Jammer  power  1  at  each  ant.  ele. 

u>j  =  Radian  frequency  of  jammer  1. 

6,  *  Uniformly  distributed  [0,2ir] 
random  phase  of  the  jammer  1 
at  the  primary  element. 

(01  Is  statistically  Independent 

of  e2.) 


> 


(Al) 


♦ 


kl 


cos(^-ak)  for  k  =  1 ,  ...  n 

I 

(2l*k-n  C0s(Vok-n)  “  2 
,  for  k  ■  n+1 ,  . . . ,  2n 


s 
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It  can  be  seen  that 


E(xjXk)  =  p1  cos(4>j1-«{»kl)  +  P2  c°s(4,j2“^k2^ 
for  j ,  k=l ,  ...»  2n 


(A2) 


The  correlation  matrix  R  is  defined  by 


R  =  CE(x.xk)] 

Substituting  from  (A1 )  and  (A2): 


(A3) 


cl+c2 

+(d1+d2> 

“(D-j  +D2) 

cl+c2 

R  = 

where  C-j ,  C2,  D-j  and  D2  are  nxn  matrices  given  by 

p.  cos j ) 

p, 

1-1.2 

j,  k  *  1 ,  ...»  n 


[Ci]jk 

£°iV 


(A4) 


Using  the  theorem  derived  in  Appendix  C,  the  eigenvalues  of  R  are 
related  to  the  eigenvalues  of  a  complex  matrix  B  given  by 

B  =  C,  +  C2  +  i[-(D1+D2)] 


(AS) 


in  the  following  manner:  n  eigenvalues  (say  ,  Xg,  ...»  XR)  of  R  are 
the  same  as  those  of  B.  The  remaining  n  eigenvalues  are  related  as 
follows: 
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*k  *  \-n  for  k  =  n+1»  •••»  2n  •  (A6) 


Writing  B  In  terms  of  vector  outer  products  of  complex  exponential 
vectors: 

B  =  p^U*  +  p2U2U*  (A7) 


where 


It  can  be  seen  that  the  rows  of  the  Hermetian  matrix  U.U.  are 

«J  J 

★ 

linearly  dependent.  Thus  UjUj  is  of  rank  1,  with  one  real  eigenvalue, 

and  all  other  eigenvalues  are  zero.  Since  the  columns  of  UjU*  are  also 

all  linearly  dependent  on  IL,  IK  spans  the  eigenspace  corresponding  to 

* 

the  one  nonzero  eigenvalue  of  UjUj. 

*  * 

Since  B  Is  a  linear  combination  of  and  U2U2»  B  be  of  ran*c 
<  2  for  n  ^  2  with  at  most  two  positive  eigenvalues.  The  eigenvectors 
of  B  will  be  In  the  space  spanned  by  Uj  and  U2-  Thus  an  eigenvector  is 
given  by 


r  -  Y,u,  ♦  y2u2 


(A9) 


where  y^  and  y2  are  complex  scalars  to  be  determined  from  the  following 
derivation. 

Combining  (A7)  and  (A9), 


BT 


>,  P,(n  V2>  *  y2U2  o2(n  *^U*U,) 


(A10) 
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Since  Br  =  Xr,  for  r  to  be  an  eigenvector 


Y2  * 

X  =  p,  n+^U, 


Y1  * 

=  <>2"+^U2Ul 


(All) 


Solving  (All),  a  quadratic  in  — ; 

Y1 


p-j  +p2  ) + 


4  (p2"P-|  )  +  P2P'|l^'j^2^ 


PlUl°2 


(A12) 


Substituting  (A12)  in  (All)  and  using  (A6), 


*1  =  ^  pi  +p2 )  + 


2 

f  (0,-P2)2  *  °1 p2 lUlU2 I 2 

m 


A 


X2  a  2^P"|+P2^  ~  4  ^1”^2^  *  Pi  P2 1  ^i  I 


n+l 


Xl;  Xn+2  =  X2 


X  =  0  for  m  Ml ,  2,  n+l,  n+2) 
m 


(A13) 


Thus  (A13)  yields  the  eigenvalues  for  the  configuration  of  Fig.  3. 

By  the  theorem  of  Appendix  C,  the  eigenvectors  of  R  can  be  expressed 
In  terms  of  the  eigenvectors  of  B,  which  will  now  be  found. 

Rearranging  (A9): 


r 


V“l  *^U2> 
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Since  the  eigenvector  can  be  linearly  scaled,  multiply  both  sides 

by  — - — -  ,  yielding  the  eigenvectors  of  the  B  matrix  of  (A5)  as 
Y1 

ri  =  C p! U1 U2 )  U1  +  blU2  (A16) 


r2  =  <P!«*U2)  U1  +  b2U2  .  (A17) 


Applying  the  theorem  of  Appendix  C,  we  can  find  the  unnormalized 
eigenvectors  of  the  R  matrix  as: 


\ 


'Re(ri)' 

r-imtr,)’ 

im(r2) 

;  Qn+1 

Redy 

corresponding  to  X-j  and  XR+1 


> 


Re(r2)' 

-Im(r2) 

o2  ■ 

Imlf2) 

II 

C\J 

+ 

C 

o 

•  A 

”R5{r27 

corresponding  to  X2  and  X^+2 


J 


(A18) 


(A18)  expresses  the  eigenvectors  for  the  configuration  of  Fig.  3-1. 
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APPENDIX  B 


The  purpose  of  this  appendix  is  to  derive  the  optimal  weight  vector, 
Wopt*  minimizes  the  output  power  of  the  antenna  array  given  in 

Fig.  1. 

The  output  power  of  this  array  (Eq.  3-3  of  the  main  text)  is: 

£  =  E[y2]  =  p  -  +  WtRW  .  (B1 ) 

Since  £  >  0,  R  must  be  positive  definite  or  positive  semi definite. 
Consequently  £  is  minimized  by  finding  the  weight  vector  which  satisfies 


IS 

3w 


0 


(B2) 


Differentiating  Eq.  B1  with  respect  to  W  gives 


=  -2P  +  2RW 


(B3) 


Therefore, 


(B4) 


Addition  of  Isotropic  white  noise  of  power  a*  to  the  system  results 


in 


V  •  t*  ♦ 'frr’p  • 


(B5) 


Since  R  Is  positive  semidefinite,  It  can  be  diagonalized  using  an 
orthonormal  set  of  eigenvectors.  Let  the  first  r  eigenvalues  be 
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non-zero,  and  denote  the  corresponding  r  unit  eigenvectors  by 
»  Qg.  ...»  Qf.  Choose  the  remaining  n-r  eigenvectors  such  that  they 
are  orthonormal,  and  span  the  null  space  of  R. 

The  matrix  which  diagonalizes  R  is 

Q  =  [Q,  Q2  ...  Qr  Qn]  •  (B6) 

Note  that 

QtQ  =  QQt  =  I  (B7 ) 

and  define 

W  =  QlW  ;  W  *  QM*  (B8) 

X'  -  QfcX  ;  X  “  QX*  .  (B9) 

Let  denote  the  white  noise  component  of  the  signal  on  element  i. 

Now, 

'V 

p  -  E[(x0  4  l^)(X  4  •  )] 

n 

-  ECx0X] 

-  E[x0QX’] 

*  QE[x0X’]  (BIO) 

P'  *=  E[x0X']  (Bll) 

P'  ■  QlP  ;  P  ■  QP'  .  (B12) 


and  define 


therefore, 
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Next,  define 


*  diagfA^,  *2*  *  *  ’  *  V*  ••**  ^ 

■  «\x« 

=  QtE[XXt]Q 

-  E[QtX  (Q^X)1] 

=  E[X* (X* )t] 

•  E[(x!)2]  . 


For  r  <  i  <  n,  we  have 


Ai  -  0  *  E(xjr]  -  0 


xj  =  0  for  r  <  1  <  n 


Hence,  the  ith  component  of  P'  is  given  by 

Pj  -  E[x0xJ] 


E[x0  0]  for  r  +  1  <  1  <  n 


»  Pj  «  0  for  r  +  1  <  1  <  n 


Combining  Eqs.  (B5,B7,B8)  gives, 
W  -  QtW 


■  «*t«„ ♦  <*** 

-  [a  +  ArV  • 


(B13) 

(BU) 


(BIS) 


(B1 6) 


(B17) 


(BIB) 


(BT9) 


From  Eq.(B18),  Pj  a  0  for  r  +  1  <  i  <  n,  so 


P^/(X1  +  oj) 


P-/(Xr  *  oj) 


If  follows  from  Eq.  (B12)  that 


and  so 


(B20) 


(B21) 


Qkp/ ( xk  +  °N}  1  5  k  <  r 


r  +  1  <  k  <  n 


(B22) 


From  Eqs.  (86, B8), 


K  *  £  Qrt  . 


(B23) 


Using  the  result  of  (B22)  gives  the  desired  expression: 


“  ■  I  ♦  «S>  *  ■ 

k«l 


(824} 


Q  E  D 
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APPENDIX  C 


Proposition 

Let  S  be  a  complex  hermitian  n  x  n  matrix  with  eigenvalues 
Aj ,  Ag,  . ...  An  and  corresponding  eigenvectors  Ej,  Eg,  . ...  ER.  Then 
a  real  symmetric  matrix  constructed  as  follows; 


R  * 


-Im(S)l 
Re  ( S  )J 


has  following  eigenvalues  and  eigenvectors 


Eigenvalues  A^ , 
Eigenvectors  k 


k 


Comments 

1.  Since  S  Is  hermitian  (1)  Aj's  are  all  real 

(11)  Ej's  all  complex 

2.  RIs  2n  x  2n  symmetric  matrix  and  as  expected  has  real  eigen 

values  and  eigenvectors. 


Proof 


R  may  alternatively  be  represented  as 


where 


R 


JLL  1  ±  £  _  |  1s_+_i7'J 
2L~ts  *  -TF !  s  +  7  J 


r_s_  i 

Jsl 

+  1 

ml 

l-1s  T 

+  2 

L-is  : 

» J 

*  J  +  J 


(Cla) 


(Clb) 

(Clc) 


(C2) 


X  L 

Consider  the  l  eigenvalue  of  S,  and  the  corresponding  eigen¬ 
vector  E^.  Construct  a  2n  dimensional  vector 


The  matrix  product 


*  *l*t  (C4) 

Thus  X^  Is  the  eigenvalues  and  4^  is  the  £***  eigenvector  of  J  for 
1  <  l ,<  n.  For  (n+1)  <  m  <  2n  we  compute  the  eigenvalues  and  eigen¬ 
vectors  of  as  follows. 

For  each  m  define  l  *  m-n.  Construct  e_  as 
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I.e. 


«»  ■  XA  ■  <«) 

thus  *m  is  the  m**1  eigenvector  with 

^  0  ~  _  • 

m  x.  m-n 

Since  4  Is  hermltlan; 

»  ^^k  ~  ^k^k  * 

Thus  the  eigenvalues  of  J  are  the  eigenvalues  of  4,  and  the  eigen¬ 
vectors  of  4  are  the  complex  conjugates  of  the  eigenvectors  of  . 

Now  construct 

Qk  "  *k+*k  (C7) 

where 

1  <  k  <  2n  . 

Using  the  definition  of  R  In  (Clc),  form  the  matrix  product  R  k* 

RQk  *  (4+*}(Sk+*k> 

*  *  rtk +  7\  ♦ 

'  Xk‘k  *  4*k  ‘‘V  Xk^k 

"Sc  ’  Vk  *  rtk  +  <“> 
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Now  consider  and  substitute  for  4  from  (C2)  and  *k  from  (C3)  or  (C5) 
depending  on  value  of  k. 

1  <  k  <  n 


(C9a) 


-ISfk-n  ♦  ,sFk.n' 

r 

;sTk-„  *  SV.  . 


•  0 

Substituting  results  of  (C9a)  and  (C9b)  In  (C8)  we  have 

“S  ■  • 

Thus  Is  the  eigenvalue  of  R  and  0^. 


(C9b) 


(CIO) 
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Eigenvalues  Ak  and  Ak+f|  are  equal,  but  the  following  argument  shows 
that  Qk  and  Qk+n  are  not  linearly  dependent. 

Consider  the  equation 


*  0  •  (Cl  3) 

If  the  only  values  (a,b)  that  satisfy  (C13)  are  (0,0),  then  ok  and 
°k+n  are  linearly  independent. 

Using  (Cll)  and  (Cl 2) 


Let 


then 


Im(&k) 


+  b 


-Im(5 


k+n-n 


Re<Wn> 


»  0 


a  -b 

Re(&k) 

b  a 

Im(6k) 

Re(Sk)  j-Im(«k) 

a 

Im(6k)  j  Re(6k) 

_  !  _ 

b 

Re(6k)  j  -Im(*kf 

_  s 

i  Re(«k) 

-  i  _ 

*  J 

det  f 
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APPENDIX  D 


Interpretation  of  the  NED  Plot  for  Small  Jammer  Separations 

As  was  discussed  in  Section  IV  of  this  report,  there  is  an  apparent 
conflict  between  the  convergence  rate  indicated  by  the  NED  plot  and  the 
convergence  rate  obtained  in  practice  when  the  angular  separation  between 
jammers  becomes  small.  The  values  obtained  from  the  NED  plot  are  large, 
showing  that  an  extremely  slow  convergence  mode  exists.  On  the  other 
hand,  intuition  and  experience  indicate  that  two  closely-spaced  jammers 
can  be  quickly  eliminated  using  a  single,  broadened  null.  This  appendix 
examines  these  two  viewpoints  more  closely  and  provides  a  resolution  of 
the  apparent  conflict. 

In  order  to  examine  convergence  rate  in  some  detail,  learning  curves 
were  plotted  for  the  three-element  triangular  array  (Figure  4-3)  using 
several  different  jammer  configurations.  Figure  D-l  is  the  learning 
curve  obtained  with  the  strong  jammer  (Power  *  100)  at  40°  and  the  weak 
jammer  (Power  *  1)  at  41°.  For  this  1°  jammer  spacing,  the  error  power 
drops  extremely  rapidly  to  a  value  of  less  than  1  x  10“^.  This  corresponds 
to  the  rapid  nulling  that  is  intuitively  expected  from  placing  a  single 
notch  on  both  signals.  Once  the  rapid  nulling  has  occurred,  however, 
the  learning  curve  becomes  essentially  flat,  and  further  reduction  of 
error  power  occurs  at  an  extremely  slow  rate.  This  is  the  slow  conver¬ 
gence  mode  Indicated  by  the  large  eigenvalue  disparity.  What  is  not 
shown  by  the  NED  plot  is  the  fact  that,  for  the  selected  initial  weights, 
the  slow  mode  becomes  Important  only  after  error  power  has  been  reduced 
to  very  low  levels. 
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Learning  Curve  for  Three-Element  Triangular  Array  with  1°  Jammer  Separation 


Learning  curves  for  the  same  array  configuration  and  jammer  powers 
are  shown  in  Figures  D-2  and  D-3  for  the  cases  involving,  respectively, 
15°  and  30°  jammer  separations.  The  15°  case  is  the  same  case  that  was 
treated  in  the  body  of  this  report.  The  NED  is  12.5,  and  convergence 
occurs  at  about  26,100  adaptations.  The  30°  case  has  an  associated  NED 
of  3.5,  and  convergence  occurs  after  some  9,800  adaptations.  The  error 
levels  that  are  reached  after  the  first  few  adaptation  cycles  differ  by 
approximately  an  order  of  magnitude,  with  the  error  in  the  15°  case 
being  the  smaller  of  the  two.  In  both  cases,  however,  error  levels  after 
the  initial  rapid  adaptation  remain  relatively  high  with  respect  to  the 
arbitrarily  selected  convergence  threshold  of  1  x  10“^.  Convergence 
along  the  slow  mode  dominates  in  both  cases,  and  the  NED  ratio  comparison 
provides  a  rough  measure  of  relative  convergence  rates: 

NED  Ratio  =  *  3.57 

Ratio  of  Adaptation  Times  =  =  2.70 

By  making  allowances  for  the  differences  in  error  levels  after  the 
first  few  adaptation  cycles,  the  NED  ratio  between  the  15°  and  30°  cases 
can  be  checked.  Crossover  of  the  two  learning  curves  occurs  at  3,340 
adaptations.  If  times  to  convergence  are  compared  from  this  point,  the 
NED  ratio  is  seen  to  accurately  reflect  comparative  rates  of  convergence 
along  the  slow  mode: 

Ratio  of  Adaptation  Times  -  g^OO  -  3340  6,460  *  3,52  * 


64 


67  x  10" 


Learning  Curve  for  the  Three- Element  Triangular  Array  with  15°  Jainner  Separation. 


It  can  be  seen  that  the  NED  provides  a  useful  measure  of 
convergence  rate  under  the  assumption  that  convergence  occurs  along 
the  slow  mode.  In  practice,  however,  convergence  will  typically  occur 
along  a  mixture  of  modes,  and  the  convergence  time  indicated  by  the  NED 
will  overbound  the  true  time.  Further  work  is  needed  to  define  the 
circumstances  under  which  the  estimate  provided  by  the  NED  is  useful  and 
to  provide  more  accurate  convergence  measures  where  the  estimate  is 
overly  pessimistic.  Consideration  is  being  given  to  an  exact  solution 
for  the  two- jammer  problem.  This  solution  will  involve  calculating  the 
optimal  weights,  determining  the  distance  from  the  initial  weights  to  the 
optimal  weights,  and  decomposing  this  distance  into  components  along  the 
eigenvectors  associated  with  the  fast  mode  and  the  slow  mode. 
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